Solar system and equivalence principle constraints 
on f{R) gravity by chameleon approach 
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We study constraints on f{R) dark energy models from solar system experiments combined with 
experiments on the violation of equivalence principle. When the mass of an equivalent scalar field 
degree of freedom is heavy in a region with high density, a spherically symmetric body has a thin- 
shell so that an effective coupling of the fifth force is suppressed through a chameleon mechanism. 
We place experimental bounds on the cosmologically viable models recently proposed in literature 
which have an asymptotic form f{R) = R — \Rc[l — (Rc/R)'^"] in the regime R ^ Re- From the 
solar-system constraints on the post-Newtonian parameter 7, we derive the bound n > 0.5, whereas 
the constraints from the violations of weak and strong equivalence principles give the bound n > 0.9. 
This allows a possibility to find the deviation from the ACDM cosmological model. For the model 
f{R) = R — XRc{R/ RcY with < p < 1 the severest constraint is found to be p < 10~^°, which 
shows that this model is hardly distinguishable from the ACDM cosmology. 
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The recent data coming from the luminosity distance 
of Supernovae la ^ , the wide galaxy surveys Q and the 
anisotropy of Cosmic Microwave Background [g] suggest 
that about 70 % of the energy density of the present uni- 
verse is composed by dark energy responsible for an accel- 
erated expansion. The cosmological constant is the most 
relevant candidate to interpret the cosmic expansion, but, 
in order to overcome its intrinsic shortcomings associ- 
ated with the energy scale, several alternative models 
such as quintessence and k-essence have been proposed 
(see Ref. U for reviews). Most of these models have the 
common feature to introduce new sources into the cos- 
mological dynamics, but, from an "economic" point of 
view, it would be preferable to develop scenarios consis- 
tent with observations without invoking extra parameters 
or components non-testable at a fundamental level. 

The simplest extension to the ACDM model is pre- 
sumably the so called f{R) gravity, where f{R) is a gen- 
eral function of the Ricci scalar R Hi (see Ref. Q for 
an early work). In Ref. Q the authors derived the con- 
ditions under which a successful sequence of radiation, 
matter and accelerated epochs can be realized. In addi- 
tion the stability conditions f^R > and /,_r_r > are 
required to avoid ghosts and tachyons for R> Ri, where 
Ri is the Ricci scalar at a de-Sitter point Q. There 
exist viable f{R) models that can satisfy both back- 
ground cosmological constraints and stability conditions 
E i, E [m, Elll, [13, E Ell • These models can satisfy 
solar system constraints under a chameleon mechanism, 
that is, a nonlinear effect arising from a lar ge d eparture 
from the background value of R [ll|, [H, Il6|| . In this 
brief report, we place constraints on viable /(i?) grav- 
ity models under the chameleon mechanism [l7[ by using 
both solar-system and equivalence principle bounds. 

We start with the following action in f{R) gravity: 



S 



-gf{R)/2 + S,n{9^.,^ra), 



(1) 



where Sm is a matter Lagrangian that depends on the 
metric and matter fields '^m- We use the unit M^^ = 
(SttG)^^ — 1, where Mpi and G are a reduced Planck 
mass and a bare gravitational constant respectively. 

We introduce a new metric variable and a scalar 
field 6, as 



(2) 



where = d f I dR. Then the action in the Einstein frame 
is given by |18| 



R/2 - (V0)V2 - V{^) 



where 



1 



V = 



Rm - f 



(3) 



(4) 



The field 4> is directly coupled to a non-relativistic matter 
with a constant coupling /3. 

In a spherically symmetric spacetime, the variation of 
the action ([3]) with respect to the scalar field cj) gives 



d^<P _^ 2 d(/) _ dFcff 
dr^ f df d(h 



(5) 



where f is the distance from the center of symmetry and 

VM<P) = V{cP)-^e^*p* . (6) 

Here p* is a conserved quantity in the Einstein frame fv^ , 
which is related with the energy density p in the Jordan 
frame via the relation p* = e^^'^p. 

We assume that a spherically symmetric body has a 
constant density p* ~ p\ inside the body (f < fc) and 
that the energy density outside the body (r > fc) is 
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p* — p*g. The mass Mc of the body and the grav- 
itational potential $c a-t the radius fj, are given by 
Mc = {4T:/3)fcP*A and <i>c = Mc/Sirrc, respectively. The 
effective potential VcS {4>) has two minima at the field val- 
ues 4iA and (l)B satisfying V^ff{<j>A) = and V^^'^ (0b) = 0, 
respectively. The former corresponds to the region with 
a high density that gives rise to a heavy mass squared 
m\ = V^g{(j)A), whereas the latter to the lower density 
region with a lighter mass squared = V^g{(j)B)- 

In the high-density regime with a heavy field mass, 
it is known that the spherically symmetric body has a 
thin-shell under the chameleon mechanism. When the 
thin-shell develops inside the body, the following thin- 
shell parameter is much smaller than the order of unity 



Arc 



(7) 



Solving Eq. ([5|) with appropriate boundary conditions, 
the field profile outside the body (f > fc) is given by 



(P[r) - — : h 0s . 
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(8) 



where the magnitude of the effective coupling, pc« = 
(3/?)(Afc/fc), is much smaller than unity when the thin- 
shell is formed. 

Let us study concrete f{R) models that can satisfy 
local gravity constraints as well as cosmological and sta- 
bility conditions. Hu and Sawicki [ll[ proposed the fol- 
lowing model 



f{R) ^R~XRc 



(R/Rc) 



2n 



(9) 



whereas Starobinsky Q proposed another viable model 



f{R) = R-XRc 1 - (1 + i?V-R,^^ 



In both models n, X and Rc arc positive constants. Since 
f{R = 0) = 0, the cosmological constant disappears in a 
fiat spacetime. Other f{R) models with similar features 
have been discussed in Refs. In these models 

a de-Sitter point responsible for the late-time acceleration 
exists at i? = i?i (> 0), where Ri is derived by solving 
the equation = 2/(i?i) Q. Note that R^ is 

not much different from the present cosmological density 
Pc ~ 10^^^ g/cm^. 

In the region R^ Rc both models ^ and pU)) behave 

as 



f{R)~R-XRc l-{Rc/R) 



(11) 



Inside and outside the spherically symmetric body 
the effective potential ^ has minima at (pA — 
-V6nXiRc/pA?''+^ and 0b ~ -V6nX{Rc/ pb?''+\ re- 
spectively. Since pa ^ Pb ^ Pc one has |0^| <C \4>b\ < 
1 and r ~ provided that n and A are not much dif- 
ferent from the order of unity. In the following we omit 



the tilde for the quantity r. From Eq. ^ the thin-shell 
parameter is approximately given by 



Arc 

Tc 



iX 
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PB 



(12) 



Let us first discuss post-Newtonian solar-system con- 
straints on the model (fTT|) . In the weak-field approxi- 
mation the spherically symmetric metric in the Jordan 
frame is 



d.s^ 



-[1 - 2A{r)]dt'^ + [1 + 2B{r)]dr^ + r^dfi^ , (13) 



where A{r) and B{r) are the functions of r. It was 
shown in Ref. [Toj that under the chameleon mechanism 
the post-Newton parameter, 7 = B{r)/A{r), is approxi- 
mately given by 



7 ■ 



1 - Arc/rc 
1 -I- Arc/vc 



(14) 



provided that the condition msr <C 1 holds on solar- 
system scales. The present tightest constraint on 7 is 
I7 - 1| < 2.3 X 10"^ [20i|, which translates into 



< 1.15 X 10" 



(15) 



For the model (fTT|) the de-Sitter point corresponds to 
A = x?"+V(2(a;f" -n - 1)), where xi = Ri/Rc- Using 
this relation together with (f>c — 2.12 x 10^^ for the Sun, 
the bound (fT51) leads to 



2n 



n 



Ri 



1) \PB 



2n+l 



— < 2.4x10" 



(16) 



For the stability of the de-Sitter point we require that 
m = Rf,RR/f,R < 1 at i? = i?i 0, which gives the con- 
(20) ditionx^" > 27i^+3?i+l. Hence the term n/2(a;f"—n—l) 
in Eq. (fTBl) is smaller than 0.25 for n > 0. Assuming that 
i?i and pb are of the orders of the present cosmological 
density 10~^^ g/cm^ and the baryonic/dark matter den- 
sity 10~^^ g/cm^ in our galaxy, respectively, we obtain 
the constraint 



n > 0.5. 



(17) 



Thus n does not need to be much larger than unity. 
Hu and Sawicki derived the Ricci scalar i? as a func- 
tion of r by considering the density profile of the Sun. 
While we have obtained the bound (fT7|) without taking 
into account such modifications, this bound is consistent 
with the one derived by Hu and Sawicki (see Eq. (67) in 
Ref. [ni). 

Let us also study the models of the type [§, [l^, [l9| 



f(R)=R-XRc{R/RcT, 0<p<l, 



(18) 



where A and Rc are positive constants. We do not 
consider the models with negative p, because they suf- 
fer from instability problems of perturbations associated 
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with negative f^RR [2i|, [2^ as well as the absence of 
the matter-dominated epoch [2^. In this case the field 
(f>B is given by (f>B = -(^/6/2)Ap(i?c//3B)^~^. Since 
the de-Sitter point, xi = Ri/Rc, satisfies the relation 
A = / (2 — p), the bound (fTS]) translates into 



P 



2 - P \PB 



R 



— \ < 4.9 X 10"" 



Taking Ri = pi = 10 g/cw? and pB 
we obtain the constraint 

p< 5 X IQ-^ 



(19) 
(20) 



Hence the deviation from the ACDM model is very small. 

Let us next place experimental bounds from a possi- 
ble violation of the equivalence principle (EP). In do- 
ing so we shall discuss the thin-shell condition around 
the Earth under the chameleon mechanism [l3|. The 
Earth has a radius = 6 x 10'^ km with a mean density 
/50 ~ 5.5 g/cm"^. The atmosphere exists in the region 
r0 < r < Tatm with a homogeneous density patm — 10"'^ 
g/cm^. The region outside the atmosphere (r > 7'atm) 
has a homogenous density pc — 10"^'* g/crn^. Defin- 
ing the gravitational potentials as <I>® = Pe'^e/S and 

^E'atm = /Oatm''atm/6, WC havC that ~ 5.5 X lO^^atm 

because p© ~ 5.5 x lO'^patm and r© ~ Tatm- Recalling the 
relation Aratm/j^atm = (0g - 0atm)/(6/3$atm), where cjiG 
and 0atm correspond to the field values at the local min- 
ima of the effective potential ([6]) in the regions r > r^tm 
and < r < Tatm respectively, we find Ar^/r© = 

-{cl)G - 0atm)/y6*e - 2.0 X 10-4(Ar-atmAatm). 

When the atmosphere has a thin-shell then the thick- 
ness of the shell (Aratm) is smaller than that of the atmo- 
sphere: rs = 10-10^ km. Taking the value ~ 10^ km 
and ratm = 6.5 x 10^ km, we obtain Aratm/'Tatm < 
1.6 X 10~^. Hence the condition for the atmosphere to 
have a thin-shell is estimated as 



Ara 



< 10" 



(21) 



Let us discuss solar system tests of EP that makes use 
of the free-fall acceleration of the Moon and the Earth 
toward the Sun. The constraint on the difference of two 
accelerations is given by 



< 10 
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a© 



(22) 



The Sun and the Moon have the thin-shells like the 
Earth [l3|, in which case the field profiles outside the 
spheres are given as in Eq. ([5]) with the replacement of 
corresponding quantities. We note that the acceleration 
induced by a fifth force with the field profile (j>{r) and 
the effective coupling (3cS is a^^*'^ = |/3off<^('')|- Then the 
accelerations a© and flMoon are [TtI 



1 + 3 



Ara 



$0 



(23) 



GMq 



O^Moon 



(24) 



where $o ~ 2.1 x 10-^ $© ~ 7.0 x 10-^° and $moo„ - 
3.1 X 10~^^ are the gravitational potentials of Sun, Earth 
and Moon, respectively. Hence the condition ((22)l trans- 
lates into 



Ara 



< 2 X lo-^ 



(25) 



which gives the same order of the upper bound as in 
the thin-shell condition (|2ip for the atmosphere. The 
constraint coming from the violation of strong equiva- 
lence principle [20| provides a bound Ar© /r© < 10"^ 
[l3| , which is weaker than . 

Let us derive constraints on the models ^ and (fTUl) 
under the bound (l25l) . On using the relation j^cl = 



V6nA(i?c/pG)'"+' > |0atm|, we obtain 



' p \ 2n+l 

PgJ 



(26) 



Taking the similar procedure we have taken to reach 
Eq. p7|) from Eq. IE]), we find the following constraint 



n > 0.9. 



(27) 



This is stronger than the bound (jl7p derived from post- 
Newtonian tests in the solar system. 

In the model (UHl) the bound leads to 

Xp{Rc/pc 



< 10 which gives the constraint 



(28) 



Thus the model is required to be very close to the ACDM 
model to satisfy the condition ([25]) . 

Let us next discuss constraints from fifth force experi- 
ments that are carried out in a vacuum pol | . Modeling a 
vacuum chamber as a sphere with radius fvaci the energy 
density is given by p{r) = for r < Tvac and p{r) = patm 
for r > 7' vac- Inside the chamber we consider two identi- 
cal bodies of uniform density pc, radius and total mass 
Mc- If these bodies have thin-shells, their field profiles 
are given by 



(r) 



47r 



(29) 



where 



is the field value when the mass squared of 



the field balances with the curvature r~^^ of the cham- 
ber. In Eq. (|29p we used the fact that the interaction 
range m^^ outside the bodies is of the order of Tvac [13] ■ 
The laboratory experiment constrains the coupling to be 
2/3^ff < 10"'^ [20|, which translates into the condition 



Arr 



< 1.7 X 10"^ 



(30) 



Thus it is crucial to have thin-shells to satisfy the exper- 
imental bound. 
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We have ^Vc/rc — —(pvac/V^^c under the condition 
that |(/)vac| is much larger than the field value \(j>A\ inside 
the bodies. A typical test body used in Hoskins et al. has 
a mass Mc ^ 40 g and a radius Tc ^ 1 cm [2^ . Hence 
the bound (l30l) translates into 



|0vac| < 10" 



-28 



(31) 



Making use of the relation A : 
Sitter point, the condition ([3T 



: x\ ^ /{2 — p) at the de- 
gives the bound 



p < 1.5 X 10" 



(36) 



For the models ^ and (fTUI) we obtain ^vac m the region 



^'■C vac 



_Qn{2n + 1)A 
Then the constraint (^1]) gives 



2n+l 
2T1 + 2 



c 



2re+l 

-l/2\ " + 1 



where C = \/&n\ 
/(2(xf 



A 



< 10 



2n + l 
2J1 + 2 



-28 



(32) 



(33) 



7r-7^5— TTTT^ ■ From the relation 

bn(2n-|-ljA Hi 

11 — 1)) we find that C is not larger 

than the order of 0.1. Using R~^''^ ~ H^^ - 10^^ cm, 
we get the following constraint: 



n > . (34) 

This is much weaker than the bounds pT]) and (|27p . 
For the model (fT8|) the field value 0vac is given by 



^ "-c ' vac 



3Ap(l-p) 



(35) 



Again this is much weaker than the bounds ([20]) and (j28p . 

In summary we have found that the models ^ and 
([To)) are consistent with the present local gravity exper- 
iments for n > 0.9, whereas the model (fT8|) is hardly 
distinguishable from the ACDM cosmology because of 
the constraint p < 10"^". These bounds are stronger 
than those derived by post-Newtonian tests in the solar 
system and arc the main results of our paper. The mod- 
els (IH) and pO|) allow the possibility to show appreciable 
deviations from the ACDM model cosmologically around 
the present epoch d, [ill, [Hj Ell • ^i^^ tie certainly of 
interest to find some signatures of modified gravity in fu- 
ture high-precision local gravity experiments such as the 
STEP [2II or GAIA [H] satellites as well as in cosmolog- 
ical observations [26| such as the galaxy power spectrum, 
Cosmic Microwave Background and weak Icnsing. 
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